Dalitz Plot Analysis of B~ — > D + iz~iz~ 



T. M. Karbach 1 , representing the BABAR Collaboration 

1 Technische Universitat Dortmund 

We present a Dalitz plot analysis of B~ — > D + ir~ir~ decays, based on a sample of about 383 mil- 
lion T(4S) — > BB decays collected by the BABAR detector at the PEP-II asymmetric-energy B Fac- 
tory at SLAC. The analysis has been published previously in pQ. We measure the inclusive branching 
fraction of the three-body decay to be B(B~ -» £> + 7r"7r") = (1.08 ± 0.03(stat) ± 0.05(syst)) x 10" 3 . 
We observe the established DJ and confirm the existence of Dq° in their decays to D + ir~ , where the 
7?2° and Dq° are the 2 + and + cm P-wave states, respectively. We measure the masses and widths 
ofDJ°andDo° to be: m D , = (2460. 4± 1.2 ± 1.2 ± 1.9) MeV/c 2 , F D * = (41.8 ± 2.5 ± 2.1 ± 2.0) MeV, 
m D * = (2297 ± 8 ± 5 ± 19) MeV/c 2 , F D > = (273 ± 12 ± 17 ± 45) MeV. The stated errors reflect the 
statistical and systematic uncertainties, and the uncertainty related to the assumed composition of 
signal events and the theoretical model. 

PACS numbers: 12.15.Hh, 11.30.Er, 13.25.Hw 

I. INTRODUCTION 

Orbitally exited states of the D meson, denoted as Dj, provide a unique opportunity to test Heavy Quark Effective 
Theory (HQET). There are expected to be four P-wave states of positive parity with the quantum numbers + (j = 
1/2), l+(j = l/2),l + (j = 3/2) and 2+(j = 3/2), which are labeled as D* , D x , D[ and D* 2 , respectively, where j is 
the sum of the spin of the light quark and the angular momentum L. Conservation of parity and angular momentum 
restricts the final states and partial waves that are allowed in the decays of the various Dj mesons. The resonances 
that decay through a D-wave are expected to be narrow (~ 30 MeV/c 2 ) and the resonances that decay through an 
S-wave are expected to be wide (a few hundred MeV/c 2 ). The D\ can only decay via a D-wave and the Dq can 
only decay via an S-wave. The D\ and D[ may decay via S-wave and D-wave. Fig. [T] shows the spectroscopy of 
the D-meson excitations and expected transitions. The Belle Collaboration has reported the first observation of the 
broad Dq and D[ mesons in B decays [3]. However, the Particle Data Group [J] considers that the J and P quantum 
numbers of the Dq and D[ states still need confirmation. 

In this analysis, we fully reconstruct the decays of B~ — > D + tt~it~ final states [2] and measure the inclusive 
branching fraction. Then we perform a Dalitz plot analysis to measure the exclusive branching fractions for B~ — > 
Dj7r~ and to study the properties of the Dj mesons. The decay B~ — > D + 7r - 7r~ is expected to be dominated by the 
intermediate states D^tt - and Dq1t~ and has a possible contribution from the B~ — > D + ir~Tr~ non-resonant (NR) 
decay. Also, the D*(2007)° (labeled as D*) may contribute as a virtual intermediate state, as well as the B* (labeled 
as B*), produced in a virtual process B~ — > B*tt~. The data used in this analysis were collected with the BABAR 
detector at the PEP-II asymmetric e + e~ storage rings. The BABAR detector is described in Q] and in the references 
therein. The sample consists of 347.23 fb^ 1 corresponding to 382.9 ± 4.2 million BB pairs taken on the peak of the 
T(AS) resonance. 

II. EVENT SELECTION 

Five charged particles are selected to reconstruct B~ — > D + ir~ir~ decays with D + — > K~tt + it + . At the T(4S") 
resonance, B mesons can be characterized by two nearly independent kinematic variables, the beam energy substituted 
mass mES and the energy difference AE: 

m ES = ^{s/2+pQ-pbY/El~p% AE = Eg — (1) 

where E and p are energy and momentum, the subscripts and B refer to the e + e _ -beam system and the B candidate 
respectively; s is the square of the center-of-mass energy and the asterisk labels the CM frame. For B~ — » D + 7r~7r~ 
signal decays, the toes distribution is well described by a Gaussian resolution function with a width of 2.6 MeV/ c 2 
centered at the B~ mass, while the AE distribution can be represented by a sum of two Gaussian functions with a 
common mean near zero and different widths with a combined RMS of 20 MeV. The AE distribution is shown in 

Continuum events are the dominant background. We suppress this background by restricting two topological 
variables: the magnitude of the cosine of the thrust angle, cos 8 t h, defined as the angle between the thrust axis of the 
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selected B candidate and the thrust axis of the rest of the event; and the ratio of the second to zeroth Fox- Wolfram 
moment [5], i?2. Small values of i?2 indicate a more spherical event shape (typical for BB events) while larger values 
indicate a 2-jet event topology (typical for qq events). We also place restrictions on m ES and AE. Then we fit the 
AE distribution to determine the fractions of signal and background events in the selected data sample. The result 
of the fit is shown in Fig. [lj it yields 3496 ± 74 signal events. 
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FIG. 1: Left: Mass spectrum [3] of cu states. The vertical bars show the widths. The dotted and dashed lines between the levels 
show the dominant pion transitions. Although it is not indicated in the figure, the two 1 + states may be mixtures of j = 1/2 
and j = 3/2, and D[ may decay via a D-wave and D\ may decay via an S-wave. Right: AE distribution. The points with 
error bars are data, the red curve is the full fit, the blue dashed curve is the background, the three shaded regions, correspond 
to the AE left side- and signal bands. The histogram shown background expected from MC. 

To distinguish signal and background in the Dalitz plot studies, we divide the candidates into three subsamples: the 
AE signal region, and two AE sidebands, all defined in Fig. [T] A background MC sample of resonant and continuum 
events is shown as the histogram in Fig. [T] There is a small amount of peaking, a fit yields 82 ± 41 peaking events. 
The background subtracted number of signal events is N S i g = 3414 ± 85, resulting in the background fraction of 
(30.4 ±1.1)%. 



III. DALITZ PLOT ANALYSIS 



In this analysis we choose the two Dtt invariant mass-squared combinations x — m 2 (D + 7T^) and y = to 2 (Z? + 7T2~) as 
the independent variables, where the two like-sign pions are randomly assigned to x and y. This has no effect on our 
analysis since the likelihood function (described below) is explicitly symmetrized with respect to interchange of the 
two identical particles. We describe the distribution of candidate events in the Dalitz plot in terms of a probability 
density function (PDF). The PDF is the sum of signal and background components and has the form: 

pnPY \ r B(x,y) [S(x,y)®1l]e(x,y) 

[ ' V) ~ /bS / D p B(x, y) dxdy + U M J up [S(x, y) ® K] e(x, y) dxdy ' [Z) 

where S(x, y) ®1Z is the signal term convolved with the signal resolution function, B{x, y) is the background term, /b g 
is the fraction of background events, and e is the reconstruction efficiency. An unbinned maximum likelihood fit to 
the Dalitz plot is performed in order to maximize the value of C — Yit^T^ PDF(xi, yi) with respect to the parameters 
used to describe S, where Xi and are the values of x and y for event i, respectively. It is difficult to find a proper 
binning at the kinematic boundaries in the cc-y-plane of the Dalitz plot. For this reason, we choose to estimate the 
goodness-of-fit x 2 in the cos 6* and m;^ in plane, which is a rectangular representation of the Dalitz plot. The helicity 
angle 9 is the angle between the momentum of the pion from the B decay and that of the pion of the Dtt system in 
the Dtt restframe; m^ in is the lesser of x and y. This analysis uses an isobar model formulation in which the signal 
decays are described by a coherent sum of a number of two-body (Dtt system + bachelor pion) amplitudes. The 
orbital angular momentum between the Dtt system and the bachelor pion is denoted as L. The total decay matrix 
element M. is then given by: 

M= PLe l ' s>I -[N L (x 7 y) + N L (y,x)} + Y,Pke l ' s>k [A k (x, y ) + A k (y,x)} 7 (3) 

L=(0,l,2) k 
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where the first term represents the S-wave (L = 0), P-wave (L = 1) and D-wave (L — 2) non-resonant contributions, 
the second term stands for the resonant contributions, the parameters p k and are the magnitudes and phases of 
the fc th resonance, while p L and $l correspond to the magnitudes and phases of the non-resonant contributions with 
angular momentum L. The functions N^x^y) and A k (x, y) are the amplitudes of non-resonant and resonant terms, 
respectively. The resonant amplitudes Ak(x,y) are expressed as A^{x,y) — Rk{m) F^ip'r') Fi(gr)Ti(p ,q, cos 9), 
where Rk{m) is the k th resonance lineshape, Fi,(p'r') and F^qr) are the Blatt-Weisskopf barrier factors [6], and 
Tl(PiQi cos 9) gives the angular distribution. The parameter m(= ,Jx) is the invariant mass of the Dtt system. The 
parameters p' , p, q and 9 are functions of x and y. The non-resonant amplitudes NL(x,y) are similar to Ak(x,y) 
but do not contain resonant mass terms. The Blatt-Weisskopf barrier factors depend on a single parameter, r' or r, 
the radius of the barrier, which we take to be 1.6 ( GeV/c) - 1, similarly to Ref. [3]. The functional forms of the Fl 
are given in Ref. pQ. For virtual D* decays, D* — > D + n~, and virtual B* production in B~ — > B*ir~, we use an 
exponential form factor in place of the Blatt-Weisskopf barrier factor, as discussed in Ref. [3]: F(z) = exp(— (z — z')), 
where z' — rp v for D* — > D + ir~ and z' = r'p v for B~ — > B*n~. Here, we set p v — 0.038 GeV/c, which gives the 
best fit. The resonance mass term Rk(m) describes the intermediate resonance. All resonances in this analysis are 
parameterized with relativistic Breit-Wigner functions: 

, / \ 2L+1 

where m and T are the values of the resonance pole mass and decay width, respectively. The terms Tl(p, q,cos9) 
describe the angular distribution of final state particles and are based on the Zemach tensor formalism [?]. The 
definitions are given in pQ. The signal function is then given by S(x,y) = \A4\ 2 . 

In this analysis, the masses of D* and B* are taken from the world averages |3] while their widths are fixed at 
0.1 MeV; the magnitude pi- and phase of the D\ amplitude are fixed to 1 and 0, respectively, while the masses 
and widths of the Dj resonances and the other magnitudes and phases are free parameters to be determined in the 
fit. The effect of varying the masses of D* and B* between 0.001 and 0.3 MeV is negligible compared to the other 
model-dependent systematic uncertainties. 

The fit fraction for the k th decay mode is defined as the integral of the resonance decay amplitudes divided by the 
coherent matrix element squared for the complete Dalitz plot: 

, _ Jpp \pk(A k (x 7 y) + A k (y,x))\ 2 dxdy 

Jk ~ J» F \M\*d X dy ■ W 

The detector has a finite resolution. For the narrow resonance D 1 ^ with the expected width of about 40 MeV, the 
signal resolution needs to be taken into account. We study the resolution on MC simulated events. We find the 
resolution to be independent of cos 9 for truth-matched events, and we describe it by a sum of two Gaussian functions 
with a common mean. The signal resolution for an invariant mass of the D~k combination around the region is 
about 3MeV/c 2 . There is a small self crossfeed (SCF) component, which varies from 0.5% to 4.0% with cos 9. For 
this component, also the resolution varies with cos(?, which we take into account. We also check the estimated biases 
in the fitted parameters due to uncertainties in the signal resolution functions are small. 

The signal term is modified in order to take into account the particle detection efficiency. Since different regions of 
the Dalitz plot correspond to different event topologies, the efficiency is not expected to be uniform over the Dalitz 
plot. We determine the efficiency function, e(x, y), by fitting twice a large sample of B~ —> D + tt~it~ MC: before and 
after the final selection was applied. The properly normalized ratio of the fit functions gives e(x,y). The efficiency is 
flat in the center of the Dalitz plot, and drops close to its boundaries. 

The background distribution is modeled using an analytic function describing MC background events. Since we 
find the Dalitz plot distributions of AE sideband events in data and in MC to be consistent within their statistics, 
we are confident that the MC simulation can accurately represent the background in the signal region. 



IV. PHYSICS RESULTS 



The total B — > D + ir ir branching fraction is calculated using the relation: B = N s i gna \/(e ■ B(D + ) ■ 2N(B + B )), 
where N aigaal = 3414 ± 85, e is the average efficiency, B{D+) = B{D+ -» I<-tt+it + ) = (9.22 ± 0.21)%, and the 
total number of B + B~ events N(B + B~) = (197.2 ± 3.1) x 10 6 . The measured total branching fraction is B{B- 
D+tt-t:-) = (1.08±0.03(stat)) x 10~ 3 . 

The Dalitz plot distribution for data is included in pQ. Since the composition of events in the Dalitz plot and 
their distributions are not known a priori, we have tried a variety of different assumptions. In particular, we test the 
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inclusion of various components, such as the virtual D* and B* as well as S-, P-, and D-wave modeling of the non- 
resonant component, in addition to the expected components of D 27 Dq, and background. We choose as the nominal 
fit model the one with the D 2 , Dq, D*, B*, and P-wave non-resonant components considered. It produces the best fit 
quality with the smallest number of components. The P-wave non-resonant component is an addition to the fit model 
used in the previous measurement from Belle |3|. The detailed fit results are: mp* = 2460.4 ± 1.2, To* = 41.8 ± 2.5, 
m Dl = 2297 ± 8, Tn* = 273 ± 12, f D . = 32.2 ± 1.3, $ D; = 0.0 fixed, f D$ = 62.8 ± 2.5, $d * ='-2.07 ± 0.06, 
f D , = 10.1 ±1.4, $ D . =3.00 ±0.12, f B : =4.6 ±2.6, $ B . =2.89 ±0.21, / P _ NR = 5.4 ± 2.4, $ P _ NR = -0.89 ± 0.18, 
/b g = 30.4 fixed, where masses are in units of MeV/c 2 , widths in MeV, fractions in %, and angles in radians. All 
errors are statistical only. The total x 2 over degrees of freedom is 220/153. The details of the other fit models in 
question are detailed in Ref. [T] . Ref. [5] argues for an addition of a Dir S-wave state near the Dtt system threshold 
to the model of the Dtttt final state. We have performed according tests [il! , which all yielded worse fit qualities than 
the nominal fit. 
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The nominal fit model results in the following branching fractions: B(B 
0.2) x 10~ 4 and B(B- D*n-) x B(D* D+tt") = (6.8 ± 0.3) x 10" 4 , where the errors are statistical only. 
Fig. [2^,-c show the rn^^Dn), m^^D-jr) , and to 2 (7T7t) projections, respectively, while Fig. [3^, and[3]D show the cos# 
distributions for the Dq and D 2 mass regions, respectively. The distributions show good agreement between the data 
and the fit. The angular distribution in the D 2 mass region is clearly visible and is consistent with the expected 
D-wave distribution of | cos 2 6 — 1/3 1 2 for a spin-2 state. In addition, the Dq signal and the reflection of D 2 can be 
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easily distinguished in the TO 2 lin (D7r) and m 2 nax (D7r) projection, respectively. The lower edge of m 
described with the D* component included than without. 
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FIG. 2: Result of the nominal fit to the data: projections on (a) mmin(Dn), (b) m^ ax (D-7r), and (c) ra 2 (7T7r). The points are 
the data, the solid curves represent the nominal fit. The shaded areas show the D\ contribution, the dashed curves show the 
Dq signal, the dash-dotted curves show the D% and B* signals, and the dotted curves show the background. 




FIG. 3: Result of the nominal fit to the data: the cos 8 distributions for (a) 4.5 < m 2 (Dir) < 5.5GeV 2 /c 4 region and (b) 
5.9 < m 2 (Div) < 6.2GeV 2 /c 4 region. The points with error bars are data, the solid curves represent the nominal fit. The 
dashed, dash-dotted and dotted curves in (a) show the fit of hypotheses 2-4 in Table [I] respectively. The shaded histograms 
show the cos 9 distributions from AE sidebands in data. 

Table [i] shows the NLL and x 2 /NDF values for the nominal fit and for the fits with the broad resonance Dq excluded 
or with the J p of the broad resonance replaced by other quantum numbers. In all cases, the NLL and % 2 /NDF values 
are significantly worse than that of the nominal fit. Fig. [3^, illustrates the helicity distributions in the Dq mass region 
from hypothesis 2-4; clearly the nominal fit gives the best description of the data. We conclude that a broad spin-0 
state Dq is required in the fit to the data. The same conclusion is obtained when performing the same tests using the 
alternative non- nominal fit models. 
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Hypothesis Model NLL xV NDF 

nominal fit 22970 220/153 

1 Dl, D* v , B*, P-NR 23761 1171/143 

2 D%,Dl,Bl, P-NR, (2+) 23699 991/144 

3 L>5, D*, B*, P-NR, (1") 23427 638/135 

4 D* 2 , D*, B*, P-NR, S-NR 23339 652/157 



TABLE I: Comparison of the models with different resonance composition. The labels, S-NR and P-NR, denote the S- and 
P-wave non-resonant contributions. 

The systematic uncertainties under consideration are detailed in Ref. pQ . The systematic effects considered include 
the number of B + B~ events, tracking efficiencies, particle identification, uncertainty on the background shapes, 
external D + branching fraction, and fit bias. 

V. SUMMARY 

In conclusion, we measure the total branching fraction of the B~ — > D + ir~ir~ decay to be B(B~ — » D + ir~Tr~) = 
(1.08±0.03(stat) ±0.05(syst)) x 10~ 3 . Analysis of the B~ — » D + tt~it~ Dalitz plot using the isobar model confirms the 
existence of a narrow D 2 and a broad Dq resonance as predicted by HQET. The mass and width of D 2 are determined 
tobemuj = (2460.4± 1.2(stat) ± 1.2(syst) ± 1.9(mod)) MeV/c 2 , Y Dl = (41.8± 2.5(stat) ± 2.1(syst) ± 2.0(mod)) MeV, 
while of the D* they are: m D * = (2297 ± 8(stat) ± 5(syst) ± 19(mod)) MeV/c 2 , T D . = (273 ± 12(stat) ± 17(syst) ± 
45 (mod)) MeV, where the third uncertainty is related to the assumed composition of signal events and the Blatt- 
Weisskopf barrier factors. The measured masses and widths of both states are consistent with the world averages [3] 
and the predictions of some theoretical models (see references in pQ). We have also obtained exclusive branching 
fractions for D* 2 and D* Q production: B(B- D\tx-) x B{D* 2 D+n-) = (3.5 ± 0.2(stat) ± 0.2(syst) ± 0.4(mod)) x 
1(T 4 , B(B- -> D*tt-) x B(D* D+tt-) = (6.8 ± 0.3(stat) ± 0.4(syst) ± 2.0(mod)) x 1(T 4 . Our results of the 
masses, widths and branching fractions are consistent with but more precise than previous measurements performed 
by Belle [S]. The relative phase of the scalar and tensor amplitude is measured to be = —2.07 ± 0.06(stat) ± 
0.09(syst) ± 0.18(mod) rad. 
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